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Abstract. In this paper we present an algorithm for computing a matrix rep- 
resentation for a surface in P 3 parametrized over a 2-dimensional toric variety 
ST. This algorithm follows the ideas of BDD09 and it was implemented in 
Macaulay2 |GS| . We showed in BDDOE^ that such a surface can be represented 
by a matrix of linear syzygies if the base points are finite in number and form 
locally a complete intersection, and in |Bot09l we generalized this to the case 
where the base locus is not necessarily a local complete intersection. The key 
point consists in exploiting the sparse structure of the parametrization, which 
allows us to obtain significantly smaller matrices than in the homogeneous 
case. 



1. Introduction 

Let J be a two-dimensional projective toric variety, and let / : ^--->P 3 be 
a generically finite rational map. Hence, S := im(/) C P 3 is a hypersurface. 
In [BDD09] and |Bot09| we showed how to compute an implicit equation for S, 
assuming that the base locus X of / is finite and locally an almost complete in- 
tersection. The work in |BDD09| and |Bot09j is a further generalization of the 
results in [BJ03, BC05, Cha06 ( BD07 on implicitization of rational hypersurfaces 
via approximation complexes. 

We showed in BDD09 how to compute a symbolic matrix of linear syzygies M, 
called representation matrix of 5, with the property that, given a point p G P 3 , the 
rank of M(p) drops if p lies in the surface S. When the base locus X is locally a 
complete intersection, we get that the rank of M(p) drops if and only if p lies in 
the surface S. 

We begin by recalling the notion of a representation matrix. 

Definition 1.1. Let S C P" be a hypersurface. A matrix M with entries in 
the polynomial ring Kpo, . . . ,T n ] is called a representation matrix of S if it is 
generically of full rank and if the rank of M evaluated in a point p of P n drops if 
and only if the point p lies on S. 

It follows immediately that a matrix M represents S if and only if the greatest 
common divisor D of all its minors of maximal size is a power of a homogeneous 
implicit equation F G K[T , . . . , T n ] of S. When the base locus is locally an almost 
complete intersection, we can construct a a matrix M such that D factors as D = 
F S G where <5 G N and G G K[T , ...,T n }. In |Bot09j Sec. 3.2], we gave a description 
of the surface (D = 0) In this paper we present an implementation of our results 
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in the computer aided software Macaulay2 |GSj .From a practical point of view our 
results are a major improvement, as it makes the method applicable for a wider 
range of parametrizations (for example, by avoiding unnecessary base points with 
bad properties) and leads to significantly smaller representation matrices. 

There are several advantages of this perspective. The method works in a very 
general setting and makes only minimal assumptions on the parametrization. In 
particular, as we have mentioned, it works well in the presence of "nice" base points. 
Unlike the method of toric resultants (cf. for example [KD06]), we do not have to 
extract a maximal minor of unknown size, since the matrices are generically of full 
rank. The monomial structure of the parametrization is exploited, in |Bot09| we 
defined 

Definition 1.2. Given a list of polynomials /o, . . . , f r , we define 

r 

Af(fo,...,f r ) :=conv(UW;)), 

i=0 

the convex hull of the union of the Newton polytopes of fi , and we will refer to this 
polytope as the Newton poly tope of the list fo, . . . , f r . When / denotes the rational 
map defining 6>, we will write Af(f) '■= Af{fi, $2, /3, A), and we will refer to it as 
the Newton polytope of /. 

In this terms, in our algorithm we fully exploit the structure of J\f(f), so one 
obtains much better results for sparse parametrizations, both in terms of compu- 
tation time and in terms of the size of the representation matrix. Moreover, it 
subsumes the known method of approximation complexes in the case of dense ho- 
mogeneous parametrizations. One important point is that representation matrices 
can be efficiently constructed by solving a linear system of relatively small size (in 
our case dimK(A„ +( i) equations in 4dimK(A„) variables). This means that their 
computation is much faster than the computation of the implicit equation and they 
are thus an interesting alternative as an implicit representation of the surface. 

On the other hand, there are a few disadvantages. Unlike with the toric resultant 
or the method of moving planes and surfaces, the matrix representations are not 
square and the matrices involved are generally bigger than with the method of 
moving planes and surfaces. It is important to remark that those disadvantages are 
inherent to the choice of the method: A square matrix built from linear syzygies 
does not exist in general and it is an automatic consequence that if one only uses 
linear syzygies to construct the matrix, it has to be bigger than a matrix which 
also uses entries of higher degree (see |BCS09| ). The choice of the method to use 
depends very much on the given parametrization and on what one needs to do with 
the matrix representation. 

2. Example 

Example 2.1. Here we give an example, where we fully exploit the structure of 
J\f(f). Take (/i, / 2 , / 3 , h) = {st 6 + 2, st 5 - 3st 3 , st 4 + 5s 2 i 6 , 2 + s 2 t 6 ). This is a 
very sparse parametrization, and we have in this case, there is no smaller lattice 
homothety of Af(f) (cf. [BDD09( lBot09j for a wider discussion on this subject). 
The coordinate ring is A — K[X , . . . ,X 5 ]/J, where J = (Xf — X 2 X t , X 2 X 3 — 
X^X^X 2 — XiX 3 ,X 2 — X X 5 ) and the new base-point-free parametrization g is 
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given by (51, <? 2 , <?3, ,g 4 ) = (2X a + X i ,-3X 1 +X 3 ,X 2 + 5X 5 ,2X + X 5 ). The Newton 
polytope looks as follows. 




For vq = 2d = 2 we can compute the matrix of the first map of the graded 
piece of degree v$ of the approximation complex of cycles (Z,) Vo (see for example 
BDD091 Sec 3.1]), which is a 17 x 34-matrix. The greatest common divisor of the 
17-minors of this matrix is the homogeneous implicit equation of the surface; it is 
of degree 6 in the variables 

Ti, . . . , T 4 : 2809Ti 2 T 2 4 + 124002T 2 6 - 5618Ti 3 T 2 2 T 3 + 668I6T1 T 2 4 T 3 + 2809T 1 4 T 3 2 

-50580Ti 2 T 2 2 r 3 2 + 86976T 2 4 r 3 2 + 212Tf Tf - 14210Ti T 2 2 T 3 3 + 3078T?T$ 
+ 13632T 2 2 T 3 4 + 116TiT 3 5 + 841T 3 6 + 14045^^^4 - 169849Ti T 2 4 T 4 
-14045r 1 4 r 3 r 4 + 261327T 2 T 2 2 T 3 T 4 - 468288T 2 4 T 3 T 4 - 7208T 1 3 T 3 2 T 4 
+157155TiT 2 2 T 3 3 r 4 - 31098^71^ - 129215T 2 2 T 3 3 T 4 - 4528Ti T 3 4 T 4 
-12673T 3 5 T 4 - 166957; 2 T 2 2 T 4 2 + 169600T 2 4 T 4 2 + 30740T?T 3 Ti 
-433384Ti T 2 2 T 3 T 4 2 + 82434T 1 2 T 3 2 T 4 2 + 269745T 2 2 T 3 2 T 4 2 + 36696Ti T 3 3 T 4 2 
+63946T 3 4 T 4 2 + 2775Ti T 2 2 T 4 3 - 19470T 1 2 T 3 T 4 4 + 177675T 2 2 T 3 T 4 3 
-85360Ti T 3 2 T 4 3 - 109490T 3 3 T 4 3 - 125T 2 2 r 4 4 + 2900Ti T 3 T 4 4 + 7325T 3 2 T 4 4 
-125T 3 T| 

In this example we could have considered the parametrization as a bihomoge- 
neous map either of bidegree (2,6) or of bidegree (1,3), i.e. we could have chosen 
the corresponding rectangles instead of M(f). This leads to a more complicated 
coordinate ring in 20 (resp. 7) variables and 160 (resp. 15) generators of J and to 
bigger matrices (of size 21 x 34 in both cases). Even more importantly, the para- 
metrizations will have a non-LCI base point and the matrices do not represent the 
implicit equation but a multiple of it (of degree 9). Instead, if we consider the map 
as a homogeneous map of degree 8, the results are even worse: For v$ = 6, the 
28 x 35-matrix M VQ represents a multiple of the implicit equation of degree 21. 

To sum up, in this example the toric version of the method of approximation 
complexes works well, whereas it fails over P 1 x P 1 and P 2 . This shows that the 
extension of the method to toric varieties really is a generalization and makes the 
method applicable to a larger class of parametrizations. 

Interestingly, we can even do better than with Af(f) by choosing a smaller poly- 
tope. The philosophy is that the choice of the optimal polytope is a compromise 
between two criteria: keep the simplicity of the polytope in order not to make the 
the ring A too complicated, and respect the sparseness of the parametrization (i.e. 
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keep the polytope close to the Newton polytope) so that no base points appear 
which are not local complete intersections. 

So let us repeat the same example with another polytope Q 7 which is small 
enough to reduce the size of the matrix but which only adds well-behaved (i.e. 
local complete intersection) base points: 




1 



The Newton polytope Af(f) is contained in 2 • Q, so the parametrization will factor 
through the toric variety associated to Q, more precisely we obtain a new parame- 
trization defined by 

(9i, 92,93, 9i) = (2Xl + W 4l -3X X i + X 2 X±, X X X± + 5X|, 2X 2 + Xj) 

over the coordinate ring A = K[Xq, . . . , X4]/ J with J = (X 2 — X1X3, X\X 2 — 
X X 3 ,Xf — X X 2 ). The optimal bound is f = 2 and in this degree the implicit 
equation is represented directly without extraneous factors by a 12 x 19-matrix, 
which is smaller than the 17 x 34 we had before. 



3. Implementation in Macaulay2 

In this section we show how to compute a matrix representation and the implicit 
equation with the method developed in [BDD09 and Bot09], using the computer 
algebra system Macaulay2 GSJ. As it is probably the most interesting case from 
a practical point of view, we restrict our computations to parametrizations of a 
toric surface. However, the method can be adapted to the n-dimensional toric case. 
Moreover, we are not claiming that our implementation is optimized for efficiency; 
anyone trying to implement the method to solve computationally involved examples 
is well-advised to give more ample consideration to this issue. For example, in the 
toric case there are better suited software systems to compute the generators of the 
toric ideal J, see |4ti| . 

First we load the package "Maximal minors^' 
il : load "maxminor .m2" 

Let us start by defining the parametrization / given by (/1, . . . , fi). 



±2 


: S=QQ[s,u 


,t,v] ; 


i3 


: el=2; 




i4 


: e2=6; 




i5 


: fl=s*u*t 


~6+2*u~2*v~6 




6 


2 6 


o5 


= s*U*t + 


2u v 


i6 


: f2=s*u*t 


~5*v-3*s*u*t ~3*v~3 




5 


3 3 


06 


= S*U*t V 


- 3s*u*t v 


i7 


: f3=s*u*t 


~4*v~2+5*s~2*t~6 



The package "maxminor. m2" for Macaulay2 can be downloaded from the webpage 
http : //mate . dm.uba. ar/~nbotbol/maxminor .m2. 
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2 6 4 2 

07 = 5s t + s*u*t v 

18 : f4=2*u~2*v~6+s~2*t~6 

2 6 2 6 

08 = s t + 2u v 

We construct the matrix asasociated to the polynomials and we relabel them in 
order to be able to automatize some procedures. 

19 : F=matrix{{f 1 ,f 2 ,f 3,f 4}} 

09 = I sut6+2u2v6 sut5v-3sut3v3 5s2t6+sut4v2 s2t6+2u2v6 I 

1 4 



o9 : 


Matrix 


ilO 


: f. 


.l=f 1 


ill 


: f. 


_2=f2 


il2 


: f. 


.3=f3 


il3 


: f. 


.4=f4 



We define the associated affinc polynomials FF_i by specializing the variables u 
and v to 1. 

114 : for i from 1 to 4 do ( 

FF_i=substitute (f _i , {u=>l , v=>l}) ; 

) 

We just change the polynomials FF_i to the new ring 52. 

115 : S2=QQ[s,t] 
ol5 = S2 

ol5 : PolynomialRing 

116 : for i from 1 to 4 do ( 

FF_i=sub(FF_i,S2) ; 

) 

The reader can experiment with the implementation simply by changing the 
definition of the polynomials and their degrees, the rest of the code being identical. 
We first set up the list st of monomials s 1 ^ of bidegree (e^, e' 2 ). In the toric case, 
this list should only contain the monomials corresponding to points in the Newton 
polytope Af'(f). 

117 : use S; 

118 : st=-Q; 

119 : for i from 1 to 4 do ( 

st=join(st .flatten entries monomials f _i) ; 

) 

120 : l=length(st)-l; 

121 : k=gcd(el,e2) 
o21 = 2 

We compute the ideal J and the quotient ring A. This is done by a Grobner 
basis computation which works well for examples of small degree, but which should 
be replaced by a matrix formula in more complicated examples. In the toric case, 
there exist specialized software systems such as |4ti| to compute the ideal J. 

124 : SX=QQ [s ,u, t , v, w,x_0 . .x_l ,MonomialOrder=>Eliminate 5] 
o24 = SX 

o24 : PolynomialRing 

125 : X={>; 

126 : st=matrix {st}-; 
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1 8 

026 : Matrix S < S 

127 : F=sub(F,SX) 

027 = I sut6+2u2v6 sut5v-3sut3v3 5s2t6+sut4v2 s2t6+2u2v6 I 

1 4 

027 : Matrix SX <— SX 

128 : st=sub(st,SX) 

028 = I sut6 u2v6 sut5v sut3v3 s2t6 sut4v2 s2t6 u2v6 I 

1 8 
o28 : Matrix SX < SX 

129 : te=l; 

130 : for i from to 1 do ( te=te*x_i ) 

131 : J=ideal(l-w*te) 

o31 = ideal (- w*x xxxxxxx +1) 
01234567 

o31 : Ideal of SX 

132 : for i from to 1 do ( 

J=J+ideal (x_i - st_(0,i)) 
) 

133 : J= selectlnSubringd ,gens gb J) 

o33 = I x_4-x_6 x_l-x_7 x_3~2-x_6x_7 x_2x_3-x_5~2 x_0x_3-x_2x_5 



x_2~2-x_0x_5 x_5~3-x_0x_6x_7 x_3x_5~2-x_2x_6x_7 I 
1 8 

033 : Matrix SX < SX 

134 : R=QQ[x_0. .x_l] 

034 = R 

034 : PolynomialRing 

135 : J=sub(J,R) 

035 = I x_4-x_6 x_l-x_7 x_3~2-x_6x_7 x_2x_3-x_5~2 x_0x_3-x_2x_5 



x_2~2-x_0x_5 x_5~3-x_0x_6x_7 x_3x_5~2-x_2x_6x_7 I 
1 8 

035 : Matrix R < R 

136 : A=R/ideal(J) 

036 = A 

036 : QuotientRing 

Next, we set up the list ST of monomials s l P of bidegree (ei,e 2 ) and the list 
X of the corresponding elements of the quotient ring A. In the toric case, this list 
should only contain the monomials corresponding to points in the Newton polytopc 

W)- 

137 : use SX 

037 = SX 

o37 : PolynomialRing 

138 : ST={>; 

139 : X={}; 

140 : for i from to 1 do ( 

ST=append(ST,st_(0,i)) ; 
X=append(X,x_i) ; 

) 

We can now define the new parametrization g by the polynomials g\, . . . ,g±. 

141 : X=matrix {X}; 
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1 8 

041 : Matrix SX <— SX 

142 : X=sub(X,SX) 

042 = I x_0 x_l x_2 x_3 x_4 x_5 x_6 x_7 I 

1 8 

042 : Matrix SX < SXX=matrix {X>; 

143 : (M , C) =coef f icients (F , Variables=>{s_SX , u_SX , t_SX , v_SX} , Monomials=>ST) 

043 = (I sut6 u2v6 sut5v sut3v3 s2t6 sut4v2 s2t6 u2v6 I , {8} I 1 I ) 

{8} I 1 
{8} I 1 0| 
{8} I -3 I 
{8} I 0| 
{8} I 10 1 
{8} I 5 11 
{8} I 2 2 1 

043 : Sequence 

144 : G=X*C 

044 = I x_0+2x_7 x_2-3x_3 x_5+5x_6 x_6+2x_7 I 

1 4 

044 : Matrix SX <— SX 

145 : G=matrix{{G_(0,0) ,G_(0,1) ,G_(0,2) ,G_(0,3)}}- 

045 = I x_0+2x_7 x_2-3x_3 x_5+5x_6 x_6+2x_7 I 

1 4 

045 : Matrix SX <— SX 

146 : G=sub(G,A) 

046 = I x_0+2x_7 x_2-3x_3 x_5+5x_6 x_6+2x_7 I 

1 4 
o46 : Matrix A < A 

In the following, we construct the matrix representation M. For simplicity, we 
compute the whole module Z\, which is not necessary as we only need the graded 
part (Zi) U(] . In complicated examples, one should compute only this graded part by 
directly solving a linear system in degree i/q. Remark that the best bound nu = z/ 
depends on the parametrization. 



i47 


: use A 


o47 


= A 


o47 


: QuotientRing 


i48 


: Z0=A~1; 


i49 


: Zl=kernel koszul(l,G) 


i50 


: Z2=kernel koszul(2,G) 


i51 


: Z3=kernel koszul(3,G) 


i52 


: nu=-l 


o52 


= -1 


i53 


: d=l 


o53 


= 1 


i54 


: hfnu = 1 


o54 


= 1 


i55 


: while hfnu != do ( 




nu=nu+l ; 



hfZOnu = hilbertFunction(nu,Z0) ; 
hfZlnu = hilbertFunction(nu+d, Zl) ; 
hfZ2nu = hilbertFunction(nu+2*d,Z2) ; 
hfZ3nu = hilbertFunction(nu+3*d,Z3) ; 
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hfnu = hfZ0im-hfZlim+hfZ2nu-hfZ3im; 

); 

156 : mi 

056 = 2 

157 : hfZOnu 

057 = 17 

158 : hfZlnu 

058 = 34 

159 : hfZ2nu 

059 = 23 

160 : hfZ3nu 

060 = 6 

161 : hfnu 

061 = 

162 : hilbertFunction(nu+d,Zl)-2*hilbertFunction(nu+2*d,Z2)+ 

3*hilbertFunct ion (nu+3*d , Z3) 

062 = 6 

163 : GG=ideal G 

063 = ideal (x + 2x , x - 3x , x + 5x , x + 2x ) 

72 35 66 7 

063 : Ideal of A 

164 : GGsat=saturate(GG, ideal (x_0..x_l)) 

064 = ideal 1 

064 : Ideal of A 

165 : degrees gens GGsat 

065 = «{0}}, «0}}> 

065 : List 

166 : H=GGsat/GG 

066 = subquotient ( I 1 I , I x_0+2x_7 x_2-3x_3 x_5+5x_6 x_6+2x_7 I ) 

1 

066 : A-module, subquotient of A 

167 : degrees gens H 

067 = «{0}}, «0}}> 

067 : List 

168 : S=A[T1,T2,T3,T4] 

068 = S 

068 : PolynomialRing 

169 : G=sub(G,S); 

1 4 

069 : Matrix S < S 

170 : Zlnu=super basis (nu+d, Zl) ; 

4 34 

070 : Matrix A < A 

171 : Tnu=matrix{{Tl ) T2,T3,T4}>*substitute(Zlnu,S) ; 

1 34 

071 : Matrix S < S 

172 : 

lll=matrix {{x_0 . . x_l>} 

072 = I x_0 x_7 x_2 x_3 x_6 x_5 x_6 x_7 I 

1 8 
o72 : Matrix A < A 

173 : lll=sub(lll,S) 
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o73 = I x_0 x_7 x_2 x_3 x_6 x_5 x_6 x_7 I 
1 8 

073 : Matrix S < S 

174 : 11={} 

074 = {} 
o74 : List 

175 : for i from to 1 do { ll=append(ll , 111_(0 , i) ) > 

176 : (m,M)=coef f icients (Tnu, Variables=>ll ,Monomials=>substitute (basis (mi, A) , S) ) ; 

177 : M; 

17 34 

077 : Matrix S < S 

The matrix M is the desired matrix representation of the surface S. 
We can continue by computing the implicit equation and verifying the result by 
substituting 

178 : T=QQ[T1,T2,T3,T4] 

078 = T 

078 : PolynomialRing 

179 : ListofTandO ={T1 ,T2 ,T3,T4> 

079 = {Tl, T2, T3, T4} 
o79 : List 

180 : for i from to 1 do { ListofTandO=append(ListofTandO,0) }; 

181 : p=map(T,S .ListofTandO) 

o81 = map(T,S,{Tl, T2, T3, T4, 0, 0, 0, 0, 0, 0, 0, 0» 

081 : RingMap T < S 

182 : N=MaxCol(p(M)) ; 

17 17 

082 : Matrix T < T 

183 : Eq=det(N); factor Eq 

We verify the result by substituting on the computed equation, the polynomials 
fi to / 4 . 

i85 :use S; Eq=sub(Eq,S) 

086 : S 

i87 : sub(Eq,{Tl=>G_(0,0) ,T2=>G_(0,1) ,T3=>G_(0,2) ,T4=>G_(0,3)}-) 

087 = 
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